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1. Introduction 

Inflation has become tlie standard paradigm to understand the initial conditions for 
structure formation and Cosmic Microwave Background (CMB) anisotropies. In the 
inflationary picture, primordial density and gravity-wave fluctuations are generated from 
quantum fluctuations and stretched to superhorizon scales during an early period of 
accelerated expansion of the universe [1]. Such fluctuations can trigger the formation 
of the CMB anisotropies taking place in a later radiation- or matter-dominated epoch. 
CMB observations [2, 3, 4, 5, 6] show that the cosmological perturbations are very small, 
of order 10~^ compared to the homogeneous background. This has justified the use of 
linear perturbation theory for a comparison with observations. There are however two 
good reasons to go beyond linear order. One is the possibility of observing some amount 
of non-Gaussianity in the CMB anisotropies, revealing fundamental details about the 
mechanism generating the primordial density perturbations. This issue is being deeply 
investigated at present, especially because of the high sensitivity to non-Gaussianity 
of the Planck satellite [7] and its successors. Such non-Gaussianities are sourced by 
self-interactions in the early universe and become visible at the level of second- or 
higher-order perturbation theory [8]. Another reason to go beyond linear order is to 
look at the loop corrections the self-interactions of any scalar field during an inflationay 
phase (and more interestingly of the curvature perturbation () can generate in the 
cosmological correlators, including the observationally interesting cases of the two- and 
the three-point correlation functions (or, in Fourier space, the power spectrum and the 
bispectrum, respectively). 

Loop corrections may lead to interesting effects which scale like the power of 
the number of e-folds between horizon exit of a given mode k and the end of 
inflation [9, 10, 11, 12, 13]. Recently there has been some renewed interest in loop 
corrections to the correlators of cosmological perturbations generated during an early 
epoch of inflation, stimulated by two papers of Weinberg [14, 15]. The reason is that 
one-loop corrections to the power spectrum of the curvature perturbation ( seem to 
show some infra-red divergences which scale like \n{kL), where L~^ is some infra-red 
comoving momentum cut-off. [16, 17, 18, 19]. |. However it has been discussed in 
Refs. [20, 21] (see also [22]) that such potentially large corrections do not appear in 
quantities that are directly observable. Indeed they concern variations of background 
quantities defined within a region of comoving size i not much larger than the present 
horizon H^^ when one considers how they change on scales L much larger than the size 
i. Such variations will be then associated to infrared corrections which however are not 
related to the local observable quantities. This means that loop corrections are under 

I By choosing the sniaUest possible value for L^^, L^^ = chH , where a^ is the scale factor at the 
beginning of inflation, the logarithm becomes In(fcL) = ln(afe/ai), which is proportional to the number 
of e-folds between the beginning of inflation and the time the mode k leaves the horizon (when k ~ UkH). 
This would have a dramatic impact for the control of cosmological perturbation theory and in particular 
for the comparison of inflationary predictions for the primordial density perturbations with observations. 
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control because it is justified to choose for tlie infrared momentum cut-off the scale £ 
not much larger than the present horizon Hq^. § 

Independently of this issue, up to now the computations of one-loop corrections to 
the power spectrum of curvature perturbations have neglected any contribution arising 
from the graviton excitations (or tensor perturbation modes) which are unavoidably 
generated during a quasi de-Sitter phase. In particular, when going beyond linear 
order, a mixing of scalar and tensor perturbations is inevitable as well. In this paper we 
compute for the first time the graviton effects at one-loop level to the power spectrum of 
the curvature perturbation ( in single-field models of slow-roll inflation. Our results show 
that such corrections are of the same order as those obtained by considering only the 
loops from scalar self-interactions (see for example [18, 19]). This further strengthens the 
point that a full self-consistent computation of one- (or higher-order) loop corrections to 
the correlation functions of primordial curvature perturbations must necessarily account 
for the graviton contribution inside the loops. Loop corrections in general appear to be 
small, neverthless an evaluation as precise as possible is required. However notice that, 
interestingly enough, we find also for the graviton-loop contributions terms which scale 
like ln(H). 

Most of the computation consists in deriving the quantum corrections due 
to gravitons to the two-point function of the inflaton field 6(f) around the time 
of horizon crossing. This is achieved using the in-in, also dubbed Schwinger- 
Keldysh, formalism [24, 25, 26]. Their contribution to the power spectrum of the 
curvature perturbation ( on superhorizon scales is then obtained by exploiting the 6N 
formalism [27, 28, 29, 30]. The SN formalism allows to obtain an expression for the 
curvature perturbation ( (in the uniform density gauge) in terms of the scalar field(s) 
fluctuations 

ax,t) = Y. Y {^^{x,t,)r, (1) 

where iV(")(t,t*) = (9/90*)"A^(t, t*) represents the derivative of the number of e-folds 
(between the final time t and some initial time t*) with respect to the initial value of 
the scalar field (in the spatially flat gauge). On superhorizon scales during inflation, 
Eq. (1) represents the perturbation in the number of e-folds of expansion between widely 
separated Hubble regions due to the initial fluctuations of the scalar field 0,^. In fact 
the time t* is left to be freely chosen, as long as it is set after horizon crossing, since at 
that time the perturbations will have settled to their classical value. One possibility is 
to take t to fall a few e-foldings after horizon crossing like for example in [31, 18, 19]; 
we will adopt this choice, motivated by the possibility of simplifying our analysis and 
calculations (we will provide details about this matter in section 3). 

The reader might be worried about the validity of the 6N expansion (1) when 
tensor perturbation modes 7jj are accounted for, since they introduce an additional 

§ In some special cases a non-linear resumination of the infra-red divergences discussed here is doable, 

see [23]. 
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degree of freedom. In fact the 6N formalism holds also in this case. The reason can 
be simply understood by considering that the spatial perturbed metric can be written 
as [32] hij oc a^(t)e^''(e'^)jj (where a(t) is the unperturbed scale factor, see also Eq. (3)) 
and, on superhorizon scales, the tensor modes can always be reabsorbed into a rescaling 
of the coordinates. || 

One-loop corrections to the power spectrum of ( are provided by the action of 6(j) 
computed up to 4th order. Both corrections due to scalar and graviton loops will be 
taken into account. The tensor fluctuations will enter the calculations up to second order 
(higher order tensor terms are unnecessary to our purpose). A possible classification 
for theories of inflation is the one that distinguishes between a canonical and non- 
canonical form for the kinetic term in the Lagrangian [33, 34]. The latter are particularly 
interesting for the study of non-Gaussianity [35, 36, 37]. We derive the solutions to the 
constraint equations for a general non-canonical Lagrangian but our loop calculations 
will be focused on the canonical case (we will cover the study of the tensor corrections in 
the non canonical case both for the power spectrum and the trispectrum of the curvature 
perturbations in a companion work [38]). 

The paper is organized as follows: in section 2 we calculate the complete action 
to fourth order in the spatially fiat gauge; in section 3 we study the graviton one-loop 
corrections to the power spectrum of the infiaton field perturbations; in section 4 we 
collect our results together with the one-loop corrections from scalar perturbation modes 
which have been already computed in Ref. [18, 19], to be combined in a final formula 
for the corrected power spectrum of the curvature perturbation; in section 4 we also 
draw our conclusions. Appendices are organized as follows: in Appendix A we face 
the problem of dealing with a lagrangian with time derivative interactions; Appendix 
B contains a detailed study of the slow-roll order of the different terms of the 4th 
order action together with the calculation of some of the diagrams with graviton loops; 
Appendix C collects the complete expressions for the one-loop two-vertex diagrams. 

2. Fourth order action for a scalar field in the spatially fiat gauge 

The action to third order with tensor modes was first calculated by Maldacena [39] 

with the ADM formalism [40]. Let us calculate the fourth order action using the same 

formalism. 

In the spatially flat gauge the perturbed metric is 

ds^ = -N^dt"^ + h,j {dx' + N'dt) {dx^ + N^dt) , (2) 

h., = a\t){e''),„ (3) 

where a{t) is the scale factor, •jij is a tensor perturbation with di'jij = 'ja = (traceless 
and divergenceless) and det(e''')jj = 1. Notice that repeated lower indices are summed 

II A rigorous proof that the presence of tensor modes docs not alter Eq. (1) can be found, e.g., in [29], 
and amounts to consider that {e'^)ij has unit determinant so that it docs not modify the local volume 
defined in terms of the local scale factor a(t, x) = a(t)e^. 
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up with a Kronecker delta, so di'-fij stands for S'^^di'jkj and 'jij = S'^^'-fij. 
The lapse and shift functions in equation (2) can be written as 



N = l + a, 
Nj = dj9 + pj , 

where Pj is divergenceless and a, 9 and j3j can be expanded perturbatively up to the 
desired order. We have exploited the gauge freedom to set two scalar and two vector 
modes to zero, thus leaving one scalar mode from N, one scalar and two vector modes 
from Nj and two tensor modes (the two independent polarizations of the graviton) from 
hij together with the inflaton field perturbation. 



We begin with the most general expression for the action of a scalar field minimally 
coupled to gravity 



S=- I dt£xVh \nR^^^ + 2NP + iV"i (EhE 



2 ' ' -'" ' ' ' '' 



E' 



(4) 



where P = P{X, 0) is a generic function of the scalar field and of its kinetic term 
X = ^g^^d^cpd'^cp, R^^^ is the curvature scalar associated with the three dimensional 



metric hij, 



Ej, 



^ {hij 



\/iNj-\/jNi] , 



E = h'^E.j. 



A dot indicates derivatives w.r.t. time t, all the spatial indices are raised and 
lowered with hij and units of Mpf = SttG = 1 are employed. In the canonical case 
P = X - V(0), with V = V{(I)) the inflaton potential. 



To 4th order we have 



i?(3) 



-jdi-faldi'Jal- 



(5) 



The next step consists in writing and solving the momentum and hamiltonian 
constraints in order to integrate out the shift and the lapse functions Ni and A^. This is 
done varying the action w.r.t. the shift and lapse functions respectively; the resulting 
equations are 



Vi 



A^-^ {e; - b)E 



AT-ip 



i^(3) + 2P - 4PxX 



X 



N~^ {E,jE'^ 



N'di^ 



di' 



E^) - 2P^xh'^di(f)dj(j) = 0. 



(6) 

(7) 
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2.1. Solving the constraint equations for a theory with non canonical kinetic terms 

The action to a given order n only requires the constraint equations to be solved up 
to order n — 2 [39, 37]. Therefore we will solve the constraints to second order in the 
metric and scalar field fluctuations. Let us employ the expansions 

a = ai + ^2, 

9 = 01 + 02. 

The momentum constraint at first order reads 

1 



2Hd^ai 



,., 2a2^'^'^' 



P 



x<, 



(8) 



where H = a /a is the Hubble parameter. Eq. (8) can be solved to get ai. Taking the 
derivative d^ of both sides of equation (8) and using the divergenceless condition for /3 
the result is 



tti 



,x 



2H 



(9) 



Using the solution found for ai, we find d'^Pij = 0, from which we can set Pi 



;i, = 0. 

Here d^ = S^Widj, which we will indicate in the rest of the paper also as didi, and from 
now on we define Pi = P2i for simplicity. 

The momentum constraint at second order is 



2Hdja2 - AHaidjai -djaid'^61 + —diaididjOi - -diaiiij 

-—d'^Pj + -iikdakj - -jikda'kj - -iikdj'jik + —^di9id'^jij 

— 2XP^xx4>otidj5(j) — P^x4>aidj5(t) 



= P,xdj64>64> + 2XP^xxdj 
+P^x<i,^S(l)dj6(j) . 



(10) 



The solutions are 



1 

AH 



d^aid'^61 — didjUididjOi 



+ -^d ^ [iijdidjai] 



^ ^P,XX Q-2 



Aam 
P 



2H 






didjdid'^-iij 



P,x 0-2 



2H 



d-'Y. 



{dj5(t)f + 509^50 



H 



d'^5(j)5(j) + dj8(j)dj5(j) — idjaidjdcj) + aid'^5 



:ii) 



where S = (9^ 



~r OjO^jyy^K 



and 
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— ^9^A = 2Hdja2 — AHaidjai -djaid'^Oi H — -diaididjOi 

2a^ a^ a^ 

- -K^iaiiij + -liikdakj - -ilikdiikj - -liikdj^ik 



2a2 



^A^^1^J - P,xdj 



2XP^xxd, 



+ 2XP^xx4>oiidj5(j) + P^x<t>(^idj5(j) - -. ^^^^.v^y^i^j. 
Let us now move to the hamiltonian constraint which provides 9i and 62 



(12) 



^a=», 



AXP 



X 



- AXP 



-Z-ai\+2P^6(P 
- 12H^ai, 



- sp.xx^ M^ - «i 



at first order and 



(13) 



AH 



d% = (-2ai) \AXPx-r + 20P,xxX'"-^ + 2XPx^5(j) + 8P,xxxX'"-Z 



+APxx6X'l 



'4o% 



AX (Px + 2XP 



XX, 



a^cj) 



.dAdM--Ad^o, 



1 ... „.2 



iiqdqdiOi + — {didj6 



-6H' + 2XPx + 4:X^P,xx 3at - 2a2 



+Aal (3X^P,xx + 2X^P,xxx) + % 



2XPx + IQX'Pxx + SX-'Pxxx 

2 



+2XP^x<j><p 
at second order. 



4XPx^ 
1 



iX^P 



XX<f> 



2^2 



a^(j) 



AX^Pxx - 2XP 



X 



+ 4 



1 1 AH 



(14) 



2.2. Reduction to the canonical case 

In the canonical case, at zeroth order in perturbation theory P = ^ — ^(0); so Px = 1 
and P^^n = —d"'V/d(j)'^ with all other derivatives of P being zero. The solutions above 
therefore reduce to 



«! 



2H 



(15) 



AH 



-d'^Oi = -2V^5(t) - 



2ai i-m^ + (0)' 



(16) 
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d~ 



d'^aid'^Oi — didjOididjOi 



(17) 



^ "-2 r.; a a „. 1 ^ g-2 



+ ^9 [iijdid^ai 



AHa^ 



didjOid'^'yij^ +^d "^ [djiikdj'yik 



AH 



^2^2 = 2ai 



^'-^i>% 



a^ 



a 



-H9idiS(j) - -^d^djOAdjO^ (18) 



a" 



hi^o^^' 



?>al - 2^2 ) (0^ - 6ij2 



—di6(j)di6(j) 



1 
Z 



■^da-fiqda-fig - 'iljilj + ^iiqdidqOl 



1 

2^ 



a^/3, 



a- 



:d'^aidjd'^ei -d^djaidmd'^ei + —d.maidmdjd'^ei 



a" 



a^ 



—djaid'^Oi -d^^djdiaidid^Oi + —d'^diaididjOi 



—didjaidid'^Oi + —dmdiaidmdidjdi + dmdj5(t)d„ 



-d' 



dj6(f)d 6(j) - dmS4>dm.djS(t) - -d [iijdiai] 



^^2 (a,-7fe,5,96^i) + ^^2 (d^^.Ae,) - ^92 (7i;ai7i/ 



2a2 

+ T<9^ iiikdakj) ■ 



(19) 



where V"^,^ = d'^V/d(j)'^ and 9 ^ is the inverse of the laplacian operator. Notice that the 
equations (15) through (19) agree with the results obtained in [41] if we set 7jj to zero. 



Sa = a^ dt(fx 



The Ath order contribution to the action is 

1 

24 
1 
a"2 

"T 



1 
3 



1 



«fa2--«2) (-6//^ + (^^ 



2V^aj50 + «! — -diS(f)diS(f) 



■[d^dje2didje, - 8^9,8^62 + d^p.didjOi) + -{<p{dje2 + /3j)djc 



+ 



■6idj5(f) 



ai 



2a2 



■7^^[lq^^a^^-iaq - -daliqdaliq) " -llljilj + iiqdidqO 



'-^liklkjdjS(j)diS(j) - ai'yijdj6(j)di6(j) + a2di 

-^djScj) {'Jijdi92 + 'JijPi + 'JijdiOi) - dk92iabdhlak - Pkiabdbjak 

1 1 

-^iahf^kdklah - ai{H'-iabdadb92 + iabdadb62 + iabdaPb) + ^iabdk1abdk02 
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+— ( - 8'yipdidjeidpdje2 - A'yipdidjOidpf^j - A-fipdpdjeidj(3i 
-dg9idajgdidj92 - dg92di'jjgdidj9i - Pgdi'jjgdidj9i 
-dg9idajgdi(3j - dg9idajgdj/3i + dg9idg^ijdidj92 

+ 2dg92dg^i,didj9, + 2dg9,dg^i^d,(3, + 2 (3 gdg^f i, d,d, 9 ,) . (20) 

Notice that the action to fourth order with gravitons was calculated in [42] neglecting 
the first-order graviton contributions for a calculation of the (tree-level) trispectrum 
in general single field models of inflation. In our paper this contribution cannot be 
neglected since we are calculating one-loop corrections to the power spectrum of the 
curvature perturbations. Moreover, as we will discuss in the next section, it turns 
out that loops of gravitons are not suppressed by the standard slow-roll parameters, 
e = —H/H"^ and 77 = (l/3)(V^0/if^), in comparison to the loops of scalar perturbations. 
From this point on, we will focus on the canonical lagrangian for our calculations. 

3. Study of one-loop tensor corrections 

The power spectrum for the curvature perturbation C, is defined by 

(Cl W4-;(t)) = {2nfPdk)5^'\k, + h) , (21) 

We will calculate it using the 6N formula 

The sums can be expanded to the desired order of loop corrections to the tree level 
power spectrum. Up to one loop we have 

+1.nWn(^) J d'q^s<p^.s<PgM,-,J. + ih - h) 

+^N^'^N^'^ I d'qdMS<Pk,5<PgMpMg+p^^k,). + {ki ^ h) 

+ (^ (^^'^y I d\d'p{5<PgMu^.gMpM,^,.^).. (23) 

where a star indicates evaluation around the time of horizon crossing. The tree-level 
power spectrum of the light scalar field at lowest order in slow-roll is given by 

(5^k,^K)* = i2nfPik)S^'\k, + h) = i27rf ^^S^'\ki + h) , (24) 
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where H^, is the Hubble parameter evaluated at horizon exit (when k = aH). The 
variance per logarithmic interval in k is given by V{k) = {k^/2TT'^)P{k).^ Notice that 
the first term on the right hand side of Eq. (23) includes both the tree level and the 
quantum one-loop contribution; the other terms represent the classical one-loop parts. 
The distinction between classical and quantum loops is understood as for example in 
[19]: quantum loops find their origin in the lagrangian interaction terms between the 
inflaton perturbations and the gravitational modes or from self-interaction of 6(j) and 
are present in the expextation values of 6(j) around the time of horizon crossing; classical 
loops are corrections coming from the expansion of ( using the 6N formula. Let us 
study the quantum loops. 

The power spectrum of Sep can be calculated in perturbation theory using the in-in 
formalism (see Refs. [24, 25, 26] or [14] which also contains a detailed review). In this 
formalism the expectation value of a field operator B(t) is given by 



{n\Q{t)\n) 







T e 



,ij[Hi{t')dt' 



Qi{t) 



T e 



-ij[Hi{t')dt' 







(25) 



where 1^2) represents the vacuum of the interacting theory, T and T are time-ordering 
and anti-time-ordering operators, the subscript / indicated the fields in the interaction 
picture, i.e. free fields and Hj is the interaction hamiltonian. The interaction picture 
has the advantage of allowing to deal with free fields only; the fields can be thus Fourier 
expanded in terms of quantum creation and annihilation operators 



(50(f,t) = / (fke'^"^ [afc50fc(t) + a!g^0fc(t) 



A 

where the commutators are the usual ones 



«fc' «fc5 



.%A'^J,V 



{2nf5^'^\k-k'), 

= {2^f5^^\k-k')8,., 



with all other commutators equal to zero, the index A runs over the two polarization 
states of the graviton and e^j are the polarization tensors. 



^ Using Eq. (24) and the 5N formula (23) at lowest order, one recovers the well known result that 
at tree level the power spectrum of the curvature perturbation C, is VQ{k) — H^/{mp-^^Te^^,), where one 
uses (TV^^^) = 47rG/e*. This equation confirms that for single-field models of inflation the curvature 
perturbation ^ is conserved on supcrhorizon scales. This is indeed true also at the non-linear level, see 
Refs. [32, 39, 43, 29]. On supcrhorizon scales also the tensor perturbation modes remain constant at 
the fully non-linear level [32]. 
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Figure 1. Diagrammatic representation of the one loop corrections to the power 
spectrum of d4> from scalar modes to leading (~ e'^) order in slow-roll. 




Figure 2. Ncxt-to-lcading {^ ^/e) order one loop corrections from scalar modes to 
the power spectrum of Scf). 



The equation of motion for the eigenfunctions S(f)k{t) can be derived in the 
approximation of de-Sitter space from the second-order action 



^2 



dri 



{Hvf 



{5<p'f - {dMY 



(26) 



(where dr] = dt/a{t) is the conformal time) and they are given by the well-known 
expression 

H 



Ukiv) 



'2k^ 



(1 + ik7]) e-'''". 



(27) 



In the same approximation, the eigenfunctions for the tensor modes 'jkiv) ^^^ given by 
ul = 2uk. 

Using the positive and negative path technique of the in-in formalism, the 
expectation value above can be recast in the form 



{n\e{t)\n) 







rfe,(t)e-*/o'^*'K(*')-^F(*')) 







(28) 



where the plus and minus signs indicate modified Feynman propagators, i.e. modified 
rules of contraction between interacting fields; schematically we have 



(r(0i02...0n))= E 
ij,lm,.. 



where the sum is taken over all of the possible sets of field contractions and 



(V)0+(V' 


(r/')0-(r/" 


(V)0+(V' 


(r/')0-(r/" 



G>{7]',r]")e{r]' -7]") + G<{r]',7]")e{7]" -r/'), 

G^V'^V"), 

G^{ri ,r] )Q{r] - r] ) + G^{r],7] )6(?7 - r/ ). 



i(0 


T 


J—oo ^ 


~iv)) 


0) 


^(0 T[e £^ rfV (i7/(V) - HYiv')) 


£jv"{Hnv")-HYiv"))] 


O), (30) 


and 


H 


7 = iff + i/f\ where H? 


and Hj indicate 
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In momentum space we have 

GUv',v") = Mv')ul{r]"), 

if (p is, for example, a neutral scalar field and Uk{r]) its time-dependent wave function. 
Remember that when computing a loop correction the external fields need to be treated 
like -|-fields. 

For our purposes, the exponentials in the expression (28) need to be expanded up 
to second order 



{Q\e{vmiL = 
+ 

X 

with 0(t) = %-;(r/)%-;(r/) 
respectively the third and fourth order parts of the interaction hamiltonian. One- 
loop corrections due to scalars in Hj were calculated in [18, 16]; their diagrammatic 
representation is given in Fig.l for the leading order and in Fig.2 for the next-to-leading 
order corrections. Loop of gravitons were ignored for simplicity in [18, 16], however 
they should be included since they are not slow-roll suppressed compared to loops of 
scalar modes. Their evaluation will constitute the primary goal of this paper. Let us 
then consider the terms in Hj that involve tensor modes. The third order action with 
gravitons has been calculated in [39] ; we will focus on the leading order term in slow-roll 
parameters and define 

iff (r/) ^ ^ J d'x^,,dMd,S<j). (31) 

The fourth order part is given by equation (88). Notice that some of the interaction 
terms involving the tensor modes in (88) appear with time derivatives. When something 
like this happens in a theory, the construction of the path integral formula requires 
additional care compared to the case where time derivatives appear only in the kinetic 
term of the lagrangian, since extra independent fields (the conjugate momenta) have to 
be taken into account. This problem will be dealt with in Appendix A. The conclusion 
of our analysis is that loops that involve conjugate momenta give zero contributions 
at one loop level to the power spectrum. It is also possible to show that in Eq. (88), 
of all the leading terms in the slow-roll expansion, only one will provide a non-zero 
contribution to the loop correction (see Appendix B for a detailed analysis), i.e. 

^(4) (^) ^ ^ J d\^,^^^^dMdM- (32) 
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Figure 3. Diagrammatic representation of the (tensor mode) corrections from ijj 
to the power spectrum of dcj). 




Figure 4. Diagrammatic representation of the (tensor mode) corrections from iJ) 
to the power spectrum of S(j). Notice that this diagram is not slow-roU suppressed 
compared to the one in Fig. 3 whereas this is not the case for scalar modes (sec Fig.l 
and Fig.2). 



Let us now begin with the one- loop one- vertex part of the diagram (given in Fig. 3) 
which we label with the subscript (IL, If); this can be written as [14], [15] 

We will study this in detail 

^drla\ri) j j^ J d\id^ q^d^q^S q^ x 

e"'^"'"%(^93)(^94)«fci(r/*)<(V)«fc2(^*)<(V)^gi(V)<(V)x 

6^^\h + (h)^'^''\k2 + gD^^'H^l + & + CO., (34) 

where the extra factor of 2 accounts for the number of equivalent diagrams obtained by 
permuting the field contractions, Uk{ri) is given by Eq. (27) and 

P,,k'y = k^y Y. ^M4 i<l) = ^k' sin^ e. (35) 

A,A' 

Integration and the use of the delta function lead to a simpler form 

{Sh-Sv*W^Svn)iiL,i.) = - ^5^'\h + fcl)^ / ^sin^^l'^ %'^^^^' -^*\l + ^krt*f 

X (1 — ikrj ) (1 -|- iqrj )(1 — iqrj ) + c.c, (36) 

This equation is exact except for the approximation of using the de Sitter space formula 
for the scale factor, aijf) = —j^, and evaluating the Hubble radius H{ri) at the time 
T]*. The reason why this is allowed is the following: the contribution to the integral 
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w.r.t. time from regions well before horizon crossing is negligible compared to the con- 
tribution due to the region around horizon crossing [39, 14, 15]; in addition to that, we 
are choosing rj to be just a few e-folds after horizon crossing, so we can assume that the 
Hubble radius (as well as any of the slow-roll parameters of the theory) will not undergo 
a big variation during this interval of time. The same approximation will be applied to 
the diagrams in the next section. 

We first solve the time integral. It is convenient to perform a change of variale like 
in [16], i.e. we set x = —krj and x* = —krj* so that 



Ht fd'q 



°° dx' k"^ 
k x"^ 



2i{x -X*) 



{Sh-,ivnS<P,^,ivl)iiL,i.) = S^'\k, + k,)^^j -^sin^e Im 

X (1 + ix*f{i - IX f {I + 4^')(i - 4^')1 • (2'') 

After integrating the imaginary part, we end up with the following result 



{'5*f,('r)«.*e(>)"))(.u.) = (38) 

= ^^'\k, + h)§i2^\[\{2> + x*^) / Y + ^(5 + Sx*^ + 2x*')jdqq , 

where the factor 4/3 comes from integrating with respect to the azimuthal angle 9 (notice 
that that the reference frame in momentum space has been chosen in such a way that 
the external wave vector k lies along the positive z axis). We now solve the momentum 
integrals. Both the logarithmic and the quadratic one exhibit ultraviolet divergences 
and the logarithmic part diverges also at very low momenta. Ultraviolet divergences 
can be treated as in flat space; the infrared logarithmic divergence is fixed introducing 
a momentum lower cutoff l~^ to be interpreted as a 'box size' [20, 21, 44, 45, 46] which 
can be fixed to be not much larger than the present horizon [20, 21]. As an example, 
consider the first integral of Eq. (38) which is convenient to split as follows 

r^ dq f^ dq f^ dq 



(39) 
q Ji~^ q Jk q 

where we have introduced an upper cutoff A. The first integral gives ln(/c£); the 
second integral can be renormalized introducing a counterterm — In(^), where k^ is a 
renormalization constant. As expected from looking at the diagrammatic representation 
in Fig. 3, there is no actual dependence from the external wave number k in the 
momentum integrals of Eq. (38). The final result for Eq. (38) can be written as 

{5<t^gM)^^kM))i^LM = ^S^'Hh + k2)^{S + x*')a , (40) 
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where a is a left over constant from renormalization. 

Let us now focus on the one-loop contribution from the 3rd order action with the 
gravitons (see Fig.4 for its diagrammatic representation) 



(5</'fc-;(r?*)(50fe:(r/*))(iL,2«) 



X 



T 



J —oo ^ 



dv"(Ht{v")-Hr{v" 



Hf 



T 



6<l>,^^iv*)S<l)^,iv*){A + B + C + D 



(41) 



where 



A= driHf dr] Hf , 

J —oo J —oo 

rv* , rv* „ 

B= dr] Hy dr] Hj , 

J —oo J —oo 

C = - dr]Hl dr] HY, 

J—oo J—oo 

D = - dr]HY dr] Hf. 

J—oo J—oo 



(42) 
(43) 
(44) 
(45) 



It is easy to check that B = A* and C = C* = D. We can write Eq. (41) as 



t:;(^*)%-^,(^*))(iL,2.) = A{-tf6^''\h + k2)k' fd'qsm^e 



(46) 



X 



dr] n* dr] ' 



{wf + wf + wf + wf) 



{Hr]'y J-oo {Ht]" 

where the factor sin^ 9 comes from contractions of the polarization tensors with external 
momenta [47] 

1 - I ^^^ I =K= sin^ e, (47) 



e^Mk'k' 



e 



V2 



qk 



V2 



and the wave fuctions Wf are 



wf{v\v") = %(^*X(V)%.(^*)^fc(V')Nii- ^(V)<i:_„-i(Wa(V)<(V') 



k-(j\' 



\k-q\' 



xQ{r] -r] ) + u%Jr] )u^^Jr] )u*(r])u(r] )e(r/ - r] )] , 



\k-q\' 



\k-q\ 



Wf{v\v") =Mfc(r/*)%(r/')M^(r/*)Mfc(r/")[M|^__^(r/')M*^__^(r/")Mg(V)M;(V') 



One-loop graviton corrections to the curvature perturbation from inflation 16 

so Wj'(?7 ^Tj ) = w9{ri ,?7 ) and is a real number and wf (r/ , 77 ) = wf*{ri ^rj). We will 
label the two contributions by A and C, so that the one-loop contribution with two 
vertices to the two point function will be broken into two parts 

Let's look in details at the two parts. 

(%-;(^*)%^,(^*))fiL,2.) = - 5^'\h + h)^, I ^4^e-^^>'^'{l + ^kv*f (49) 

/K J q \k — q\ 

X r ^e''''^''-''-\^-'^\l-tk7]'){l + tq7]'){l+t\k-q\7]') 
J-00 rj' 



X 

' -00 7] 



rl " 

dr^^i(fc+g+|fe_,-i)^"^^ - ikr]"){l - iqv"){l - i\k - q]v") + c. 



The second time integral has e*^'' — \r + -rj — i (^-r^j as its primitive function, where 

g = k + q+\k — q\,h= —qk — {q + k)\k — q\ and c = qk\k — q\. This should be evaluated 
between —00 and 77 . It is soon evident that the lower bound represents a problem for 
this evaluation. We need to remind ourself, though, that the choice of the integration 
time contour needs to be deformed and to cross the complex plane to account for the 
right choice of the vacuum [39]. This is done by integrating in a slightly imaginary 
direction, i.e. taking rj ^ rj + ie\r] |, where e is a fixed small real number; so for 
example 

, ikrj 

dr^e''"^ = -—. (50) 

3 ik 

With this contour prescription, our integral in 77 vanishes at —00. Performing the 
same change of variables as in (37) 

(k c I (gh — c\\ ( d I s ,2 c ri" 



v2 



(51) 



where (i = q—k+\k—q\ and s = kq+{k—q)\k—q\. The result of the integration w.r.t. time 
is a polynomial function of sin2x*, cos2x*, Si(2a;*), Ci(2x*) and their products with 
coefficients which depend on g^ b, c, d, s, k. Notice that in the large scale limit x* — > a 
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singulariry similar to the one found in [16] shows up in our result. However we evaluate 
the power spectrum of 6(f) just a few e- folds after horizon crossing, so we are safe from 
these kind of singular behaviour [19]. 

The next step consist in performing the momentum integral. The integrals we need to 
evaluate are of the following kind 

d^q sin^g 

^—-3/(^1' (52) 

y \q — k\ 

where /(g) is a sum of functions of momentum. Let us begin for simplicity by considering 
the constant term of the sum, i.e. let us study 



l' \q-kf' 



(53) 



For the specific case of equation (53) the integrand function has singularities at 
g = and ai q = k and shows no ultraviolet singularities. Based on an approximate eval- 
uation performed considering a sphere of radius i~^ around q = 0, where i~^ <^ k, the 
integral is proportional to a function ln{ki). The same result can be obtained working in 
a small sphere around q = k after a change of variables go = q—k. The contribution from 
large values of q is negligible w.r.t. the ones from the singular points, so the integral over 
the whole momentum space is expected to be proportional to ln{ki). The exact value of 
the integral can been found after a change of variable from the (g, 6) to the (g,p) space, 
where p = |g - k\ and is equal to (167r/225P) (1 + 301n(H)) ~ k'^ (10"^ + lOln(H)). 

Integrating Eq. (51) we find ultaviolet power law and logaritmic singularities in 
addition to infrared logaritmic contributions. The final result of the integration is a 
function of x* = e~^*, where N^ = ln(a*/a/c) is the number of e-foldings from horizon 
crossing 

(50fc;(^*)%-;(^*))fiL,2.) = ^S^'Hh + h)^{a, \n{k) + a, \n{ki) + a^) , (54) 

where Oi, 02 and 03 are functions of x* (see Appendix C). We are calculating the two 
point function for the scalar field a few e-foldings after horizon crossing, so x* may 
be chosen to range between 10~^ and 10~^. In this range ai ~ 0{1) and negative, 
a^ = -16/(15x*2) + (8/15)(5 - 8Ci(2x*)) and 03 = -8/(225a;2) + C(l) + p, where p is 
a left-over scheme- dependent renormalization constant of the kind present in equation 
(40). 

Let us now move to part C of Eq. (41) which we give below 

zAt J g \k — q\ 
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-n* drj 



-oo Tj 



''I* drj' 



-wv 



V 



Q*iv" 



Ht fd^q sin^e 



S''\h + h)^J^-;^^^{^ + Wr 



X 



Re / drj e 



^gr, Qijl 



'x\ 2 



Im I df] e 



tg^Qi;n} 



where Qij] ) = 1 + igrj +bq — icrj . 
Let us integrate over conformal time 

dx' P 



(55) 



X. k x'^ 



e-^^' 



g , h ,2 c /3 

1 + 1-X+—TX -l-r^X 

k k'^ k-^ 



e *fex* 
A;2 



k'^ ck igh — c)k'^ 

\ X +i 

X g g^ 



(56) 



where again the integration has been performed by continuing rj to the complex plane, 
i.e. [t] —* f] + ie\ri |), and then taking the limit e — ;> 0. 
We are now ready to integrate over momentum 



(%q(^*)%-;(^*)) 



c 

{lL,2v) 



s"Hh.Lr,)^4 



Ht fd^q sin^e 



^' \k-qf 



*2 



X 



k^ 2ec 2k%c ^c^ Pb'^ 



X 



*2 



+ 



+ 



+ 



X*^C^' 



(57) 

9 g-' 9^ g^ g^ J 

Similarly to what we have done in part A, one can check that there are no ultraviolet 
singularities in the remaining five integrals although some infrared logarithmic 
contributions are still present and the final result is 



H^ 



{S^kSv*W,^,ivl)fiLM = ^S'^'Kh + h)^ (ci + C2 ln(H)) 



(5^ 



where ci = (1/225) (8/x*^ + 107 + 50x*2) and ca = (16/15x*2) + (4/15). Notice that 
the (x*)~^ coefficients in ci and C2 exactly cancels the (a;*)~^ coefficients in 02 and 03. 
This is not surprising: based on [16, 19], we expect we might observe a logarithmic 
singularity if we push x* ^ in our results (which is indeed present in the Ci(2x*) term 
of 02), but no power-law singularities are actually expected. 

4. Final results and conclusions 



Let us now collect our results in the final formula for the power spectrum of the curvature 
perturbation ( computed up to one-loop level. This can be derived from Eq. (23), which 
follows from the SN formula, and the final expression reads [31, 19] 

{CkSt)Ck,it)) = {'27lf6^'\h + h) [ {N^'^y (Ptrocih) + Po„e-loop(A;i) 

+ N^'^N('^J^^Bik„q,\h-ql) 
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+ N^'^N('^P,Uk) J 70j^trcc(g)] , (59) 

where Ptrcc(^) is the tree level power spectrum (24) and B{ki, k2, k^) is the bispectrum 
of the scalar field defined by 

{6<I),-^6<I),-J<I)^^) = {27ry'6'^'\h + h + h)B{k,, k^, k^) . (60) 

In previous works the computations of one-loop corrections accounted just for the 
contributions coming from the infiaton field peturbations (see, for example, Refs. [16, 
17, 18, 19]. The focus of this paper is on including also the one-loop terms arising 
from interactions between the tensor (graviton) modes and the scalar field in the case 
of single-field models of slow-roll infiation. At this order the graviton contribution enter 
the final formula as a correction to the power spectrum of the infiaton field around the 
time of horizon crossing 

-< one-loop ("^j = -< scalar!, "^j + -ftcnsorl"^) ; (61) 

and as such it is due to genuine quantum effects, while it does not affect the remaining 
terms of Eq. (59), which, according to Ref. [19], can be considered as classical 
contributions arising after the perturbation modes leave the horizon. + 

Summing the main results of the previous section, Eqs. (40), (54) and (58), we 
finally find the one-loop graviton correction to the infiaton power spectrum 

{S<PkM*)^<Pk,iv*))iL = 7^S^'\h + h)^ IfiHk) + /2ln(H) + /s] , (62) 

where 

fi = -^{5 + 5x*^ + 2x*^) , (63) 

/2 = «2 + C2 , (64) 

and /a is given by a left-over scheme- dependent renormalization constant plus 
contributions of order 0{1) (see Appendix C for the complete expressions of 02, C2 and 
/s). If we calculate the two point funtion of S(f) a few e-foldings after horizon crossing, 
i.e. X* ranges for example between 10~^ and 10~^, /i reduces to a negative constant of 
order 0(1) and /2 ~ 4 (1 - Ci(2x*)) ~ 0(10). In the limit where x* -^ 1 both /i and 
/2 turn out to be of order unity. 

In order to understand which is the dominant contribution in Eq. (59) and how 
big it is, one needs to (i) know the slow-roll order of the coefficients N^^^: N^'^^ ~ e~^^^, 
iV(^) ~ e°, A^''^-' ~ e^/^; (ii) compute the integrals involving the power spetrum P{q). 
This is discussed in details in Ref. [19] (for the case of scalar perturbations only), see 

+ Notice however that the computation of (C^r (^)Cfc' (t)) could be carried out without using the 6N 
forniahsm, working in the comoving gauge and making use only of the Schwinger-Keldysh path integral 
formula as in [39]. In this respect, the two types of terms, "classical' and "genuine quantum effects", 
arc not really supposed to be entirely different. 
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in particular Sec IV of [19]. It turns out that the crucial quantity is represented by 
the number of e-foldings of inflation between the times of horizon exit of the mode i~^, 
which corresponds to the infrared cutoff, and the time of horizon exit of the mode k we 
want to observe. However, to deal with observable quantitites one has to choose i not 
much bigger than the present cosmological horizon Hq^ [20, 21]. 

The relevant point about Eq. (62) is that it gives in Eq. (59) a contribution which is 
of the same order of magnitude as those coming from loops which accounts for scalar 
perturbations only. Since in terms of the slow-roll parameters f A^*^^M ~ e*~^ the 
magnitude of the one-loop graviton correction turns out to be 

APl'^°^ik)r^^a{k)^V!{k), (65) 

where we have used Eq. (24) for the power spectrum of the infiaton field. In Eq. (65) a{k) 
includes the various coefficients of Eq. (62), and it is 0{1). Eq. (65) allows a more direct 
comparison with the results of Ref. [19], showing that the graviton contributions to the 
one-loop corrections are comparable to the ones computed only from scalar interactions. 
Notice that also for the tensor contributions we find terms of the form \n{ki). 

One-loop corrections to the power spectrum of density perturbations from infiation 
are small, still a precise and self-consistent computation requires to account also for 
the loops which are switched on by the tensor (graviton) modes. First, going beyond 
linear order, the tensor perturbation modes produced during infiation unavoidably mix 
with scalar modes. This fact alone would require to include the tensor modes for a self- 
consistent computation. Most importantly, despite a naive expectation suggested by 
the fact that the power spectrum of the tensor modes is suppressed on large scales with 
respect to that of the curvature (scalar) perturbations, our results show explicitly that 
their inclusion is necessary since their contribution is not at all negligible with respect 
to the loop corrections arising from interactions involving the infiaton field only. 
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Appendix A. Path integral formalism for a lagrangian with time-dependent 
interactions 

The propagator of two fields 0i and 02 is defined by (see for example [48]) 

(0102) = / D(j)DUe' / M^^-H) ^ (66) 
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where 11 is the momentum conjugate to and H is the hamihonian density. If H is 
quadratic in 11, as it happens for example in flat space-time for a field governed by a 
lagrangian L = J d^x (^d^d^cp — ^(0)) the square in the exponent can be completed 
and the integral in 11 evaluated and all is left is 

(0i02) = //^0e*''. (67) 

So, if interaction with time derivatives appear in the lagrangian, 11 and are 
independent fields in the path integral. This will provide some extra vertices that need 
to be accounted for in the Feynman diagrams. We sketch a derivation of these extra 
vertices. It will turn out to be similar to what Seery does in [18], although complicated 
by the presence of gravitons. To keep the calculations easier we will at first ignore spatial 
derivatives and tensor indices, this will also make the notation simpler. Also, we won't 
keep track of all of the numerical real coefficients since this is not crucial for getting to 
the final result; it is instead very important to keep track of imaginary coefficients, time 
derivatives and powers of the scale factor a factors, and we will make sure they are all 
accounted for in our analysis. 
The total action is S = J drj {L^ + L^), where 

L^ = a^7 + r^7 6(f) + r^7 + T^6(j) 7 + T^^6(j) 7 + r^^7 

+ ^0^^(50 + A(^(^^7 + A^^^7 , 
L^ = a'^S(j)'^ + Ti6(j)'^ + 12^0'^ + cj50'^ + X6(f)\ (68) 



where / = df /dr] and where we define 

r^ ~ a7, (69) 

r<^ ~ a6(P, (70) 

Ti ~ 00^50, (71) 

T2 ~ a'50^ (72) 

r^0 ~ a^507, (73) 

r^^ ~ 0^72, (74) 

u ~ a6(f), (75) 

Xcj,<jy^ ~ a50^, (76) 

A^^^ ~ a507^, (77) 

Xcj,<jry ~ 07^0^, (78) 

A^^^ ~ 07^. (79) 

Notice that in equations (69) through (79) we use the equivalence symbol meaning 
that we skeep details about integrations in momenta and real coefficients. 
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The conjugate momenta are 

n^ = -— ^ = a^7 + X-y^^ + X^^ + r^0(50 + T-yScj) + 1^7 Sep + r^^7 + r^7 , (80) 

07 



n<^ = 



5L 



/2 



5{5(P 



r^7 50 + r^07 + X^^y + a (50 + ri50 + 12^0 + A^^^ + cc;50 . (81) 



We solve perturbatively the equations (80) and (81) in order to derive 7 and 50 to 
fourth order 



7 



n^ + A^^^ + A 



777 



\(^07 



a"'r^^n^ 



a-^r^n'^n'^ + a-T^mm 






+ a-^r^<^rin^ + a-^TiX^^^ + a-'^riiirin'^ + a-^riran"^ 



12) 



(83) 



The next steps are: derive the hamihonian H = IP'j + n'^50 — L (7, 50, 7 , 50 
where we need to plug in the solution (82) and (83) for 7 and 50 ; construct the action 
as S* = 5*0 + Sn, where Sq = J drj (L^ + L^) and S'n includes the terms that depend on 
the conjugate momenta of the fields (a change of variables similar to the one that Seery 
performs in [18] over the conjugate momenta will also be necessary). 
Let's consider the vertices in Su that are involved in the corrections to the one loop 
point function for the scalar field 



S'n D drj a '^ri50 11' 



a-^riH'^n'^ + a-^ri7'n'^ + a-^isn'^n"^ 



+ a" 



n^n^ + a-''cj50 n'^n'^ 



^4) 



The first three vertices belong to the third order part of the action; a~^ri50 11'^ and 
a~^ri7 11'^ provide a correction to the two point function at one loop with two vertices. 
Because of the presence of Fi which involves a factor of 0, it is subleading in slow roll 
order w.r.t. the corrections coming from fourth order vertices. We will therefore neglect 
these diagrams. The same applies to the second vertex, a~^rin'^n'''', although this may 
contribute to correcting the one point function 



(50,^7/*)) D a 



Hi 



d^qfM f dr]'6{0) (1 - ikT]') e^^''' + c.c, (85) 



where Ci is a numerical real coefficient and 5(0) is the Dirac delta function deriving 
from the propagator of the II'^'s and /i is a scalar function of the internal momentum. 
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The main contribution to the integral is due to times around horizon crossing since at 
early times the rotation to imaginary plane of the contour integral makes the exponent 
decrease rapidly to zero and moreover 77* was chosen to be just a few e-folding after 
horizon crossing. Also, since the integrand function goes to zero as r] approaches zero, 
we get a good approximation of this integral taking the upper limit 77* — > 0. The result 
is purely imaginary and it cancels out with its complex conjugate. 
Let's now move to the fourth order vertices (Fig. 6). From a~^r2n'^n'^ we have 

i^^k, ivlS4>,-^ ivD^^C^J d'qh {q) f_^ dri'S{0) (l - ^kv) ' e''^^' +c.c. (86) 

The same consideration as in (85) apply to the integral above, which gives a zero 
contribution, as well as the following diagrams (corresponding to the last two vertices 
in (84)) 

El 

^ k^ 



^ Ci f d^qfM r dr]'6{0)r]'^ (l - ikT]') e^^'^"' + c.c. (87) 



Appendix B. Study of leading slow roll order vertices in the fourth order 
action 

We are interested in computing the correlators just a few e- foldings after the scales we 
consider cross the horizon, so we can assume that the slow roll parameters remain small 
and can be treated as constants during this length of time. It is then correct to limit 
our interest to the leading order slow-roll contribution to the action . 

Let's start from the study of the slow-roll order of the fluctuations derived as 
solution to the constraint equations: 



62 = sRg, [6<P^] + V~eSg, [50, 7] + Tg, [502] + 52 U^^ [^^2] ^ y^^ [^2] 

+ e^'^Wg,[5<P,ll 
(5, = eR,[5<f] + v^5,[50,7] + T,[5<P^] + V,Yl\ 

where 5'[50, 7] is a linear function of 50 and/or its derivatives and a linear function of 
7 and/or its derivatives, R[5(j)^] is a quadratic function of 50 and/or its derivatives and 
so on. 
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Notice that the first order fluctuations are subleading (~ ^/e) w.r.t. the second 
order ones (~ e^). This criterium allows a suppression of a large number of terms in 
the Ath order action based on keeping the leading order (i.e. ~ e^) terms only 



S, 



dtd X 



-^{d.^j + dj(5,)di(3, - ^(50 {8^62 + (3j) d,5<f) + 3H 



2^2 

On 



2 -jliklkjdjScpdi' 



a2di 



dk02iabdblak + ^iabdklabdA 



+ -^ [Pkiabdblak - -iabf^kdklab 



It can be easily shown that the terms in the action that do not contain the gravitons 
reproduce the ones in equation (37) of [41]. The contribution to the power spectrum 
due to these vertices has been calculated by these authors, but only for the scalar part. 
We will then focus on all the tensor contributions from these and from the remaining 
terms. Interaction vertices with both two and four tensor fluctuations will be obtained 
once the expressions for 02 ? (^2 and I3j are plugged in the action. The terms in the action 
that we need for constructing Feynman diagrams with one loop of gravitons are 



5*^2 



a^ / (Txdt 



1 



1 



- ^/^.5^/5. - -28H 



J ^2 - ■j-^djS4>diS(j)'yik'ykj 



-6^dj6(j)f3j + -r-^{2iabda-fak {dk62 + Pk) - iabdklab{(3k + dkO^, 

Let's plug the expressions for f]j and 62 into (88) considering the terms with two 
gravitons. The result is an ensemble of vertices which can in principle contribute 
to the one loop corrections to the power spectrum of the scalar field. Apart from 
djS(j)diS(j)'jik'jkj, all of the other terms contain time derivatives of one, two or three of 
the four fields 



- dmS^dmdj6(j)j iiikdakj - liida'kj) 
+ d~^ (dmdj6<pdmS(j) — d'^64>dj6(j) + djC 

- dm6(i)dmdj6(j))d~^ {iikd^kj - ludiikj) 



^9) 



54)3^8.92 D 



2 
2 



d "^dj ^r-^da-fiqda'Jig + iljilj 



S4>djS(f>Pj D S^djS(P—d-^dj 



iikdi'jkj - likdiikj 



^' --2o r ^0-^-2^- - '2 1 



iabdalbkdk02 D iabdalbk-TTjd Ok - QH 8 S - 50 ^5^509, 

AH L a^ 



(90) 
(91) 
(92) 
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iabdkJabPk 3 iabdklabd ^ 

- dmS4>dmdkS( 

iabdk1abdk02 D iabdklab 



dmdkUCpUm 



-QHd^Hd^ 



d' 



+ 9fc509^ 



9, 



•2 1 

6 7^di5(t)dii 

a^ 



25 



(93) 



(94) 



(95) 



We will now prove that the vertices that include time derivatives do not actually 
contribute to the two point function. First of all notice that the tensor fields carry 
polarization tensors e^j with the property q^eij = and are always contracted with 
other tensor fields in the calculations; this implies that, if a partial derivative index is 
contracted with a tensor index, that diagram will be zero. Based on this observation, 
we can ignore several of the vertices with time derivatives. We are eventually left with 
only two of them, that we will call Vi , V2 and V3 



Vi ~ dj {5(l)dj6(t)^ d ^ {da%cdaibc) , 

V2 ~ dj (^6(l)dj6(f)j d""^ iiabiab) , 



(96) 
(97) 



V3r^iabdk7ab{Pk + dke2). (98) 

where (98) is given by the sum of (94) and (95). Notice that the 'jij fields need to be 
contracted between each other and that J2\\' ^iq^iq =constant [47]; the derivatives of 
6(j) contract with derivatives of 7, so this produces k ■ q factors. Therefore we have 



(%q(^*)%-;(^*))vi+v^. ~ ^S^'\k, + g)Ht J ^h{q')k ■ qj'^ driUv 



CO., 



(99) 



where fi{q^) and /2(?7 ) are some functions of q^ and 77 . This contribution is evidently 
zero for symmetry reasons. 



Appendix C. Complete expressions of one-loop two- vertex diagrams to 
leading order 

In the following we provide the explicit expression for Eqs. (54) and (58). Equation (54) 
reads as 



Ht 



{S4>^^Sv*)S^kM))tiLM = TrS^'\h + k2)^(ai Hk) + a2Hki) + a-s 



k^ 



(100) 
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where 
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ai= -—(5 + 5x*^ + 2x*^., 
15 ^ 



(101) 



a2 



*2 



2 + 5 - SccCTc + ^iras - SasOsjx*^ - Siirac + 2aca. 



*2 



15x 

2acd-c)x*'^ + (1 + Sacac — 4:nas + 8asd's)x 



*4 



(102) 
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3120 + 15136acac - 4507rVcac - 75687ra, 



1800x*2 
+ 2257rV, + 15136(T,<5-,)x*^ - (l51367rac - 4507rVc + 30272^^(7, 

- gOOvr^acCr, - 30272o-,a, + 9007rVca,)x*^ - ( - 672 - IdlSGa^ac 

+ 4507rVeac + 75687ra, - 2257rV, - 15136(T,5-e + 4507rV,a,)x*^ 



208a; 



*6 



P, 



(103) 



and p is a constant left over from renormalization of ultraviolet divergences. We have 
defined 



o"s = sin2x*, 
(Tc = cos2x*, 
(Ts = Si(2x*), 
a, = Ci(2x*), 

where Si and Ci stand for the sine-integral and the cosine-integral functions, i.e. 



Si(x) 



^ sin(t) 



t 
^ cos(t) 



dt, 



-dt + \n{x) + 7, 



Ci(x) 

JO t 

with 7 indicating the Euler Gamma function. 
The expression for Eq. (58) is 



where 



H^ 



^^iv*)S<p>^M))'iiLM = ^5^'\h + h)-^ {c, + c^Hki)) , (104) 



k^ 



Cl 



225 Vx*2 



^r + 107 + 50x 



*2 



C2 = TT 



4/4 
15 Vx*2 



(105) 
(106) 
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Finally the quantity /s appearing in Eq. (62) is given by 

/3 = a3 + ci + a', (107) 

where a' = 2 (1 + x*^/3) a from Eq. (40). 
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